基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

自由基浓度测定方法的再验证及其在煤化学中应用

对前人建立的标准曲线法测煤中自由基浓度进行优化,以DPPH标准样品和基准样品的二次积分面积比值为新参数,结果显示新参数标准曲线法的实测值与理论值相对误差都在5%以内;重复性、复现性实验的相对标准偏差都小于3%。将新参数标准曲线法用于分析不同煤化程度煤和新疆黑山煤(HS)沥青质的自由基浓度,发现随着煤化程度增加,其煤中自由基浓度逐渐增大,从低阶褐煤的8.531×10~(17)/g上升到高阶无烟煤3.37899×10~(19)/g;而在HS煤液化过程中,随着加氢液化温度的升高,其沥青质自由基浓度逐渐下降,从290℃的1.5793×10~(18)/g降到450℃的7.410×10~(17)/g,沥青质自由基浓度变化趋势与其产率变化趋势相一致。     

A Theoretical study on the charge and discharge states of Na-ion battery cathode material,Na1+xFePO4F

Na₂FePO₄F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na₂FePO₄F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na_1.5FePO₄F has the P2₁/b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P2₁/b11 Na_1.5FePO₄F and Pbcn Na₂FePO₄F or Na₁FePO₄F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na₂FePO₄F has diffusion paths in the a- and c-axis directions, Na₂FePO₄F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.     

Shape morphing and topology optimization of fluid channels by explicit boundary tracking

An integrated shape morphing and topology optimization approach based on the deformable simplicial complex methodology is developed to address Stokes and Navier‐Stokes flow problems. The optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. The shape evolution is realized by deforming the curves while maintaining a high‐quality mesh through adaption of the mesh near the structural boundary, rather than performing global remeshing. Topological changes are allowed through hole merging or splitting of islands. The finite element discretization used provides smooth and stable optimized boundaries for simple energy dissipation objectives. However, for more advanced problems, boundary oscillations are observed due to conflicts between the objective function and the minimum length scale imposed by the meshing algorithm. A surface regularization scheme is introduced to circumvent this issue, which is specifically tailored for the deformable simplicial complex approach. In contrast to other filter‐based regularization techniques, the scheme does not introduce additional control variables, and at the same time, it is based on a rigorous sensitivity analysis. Several numerical examples are presented to demonstrate the applicability of the approach.     

The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

Local discontinuous Galerkin method for a nonlocal viscous conservation laws

The purpose of this article is to investigate high‐order numerical approximations of scalar conservation laws with nonlocal viscous term. The viscous term is given in the form of convolution in space variable. With the help of the characteristic of viscous term, we design a semidiscrete local discontinuous Galerkin (LDG) method to solve the nonlocal model. We prove stability and convergence of semidiscrete LDG method in L² norm. The theoretical analysis reveals that the present numerical scheme is stable with optimal convergence order for the linear case, and it is stable with sub‐optimal convergence order for nonlinear case. To demonstrate the validity and accuracy of our scheme, we test the Burgers equation with two typical nonlocal fractional viscous terms. The numerical results show the convergence order accuracy in space for both linear and nonlinear cases. Some numerical simulations are provided to show the robustness and effectiveness of the present numerical scheme.     

基于PI法的门源MS6.4地震前震中附近地震热点图像异常变化研究

基于1980年以来青海省和甘肃省区域地震台网的定位地震目录，分析了研究区最小完备震级，选取空间网格尺度为0.2°×0.2°、地震活动异常学习时间段和预测时间段均为3 a，应用图像信息（PI）法，分析了门源M_S6.4地震发生前震中附近地震热点图像的异常变化过程。结果表明，在2013年1月1日至2016年1月1日的预测时间窗内，冷龙岭断裂和祁连山北缘断裂附近存在明显的地震热点，且地震热点图像颜色偏深，发震概率较高，向前滑动至2014年1月1日至2017年1月1日，2015年1月1日至2018年1月1日和2016年1月1日至2019年1月1日3个预测时间窗，此两断裂附近仍存在地震热点，地震热点图像颜色逐渐由深变浅，而门源M_S6.4地震就落在PI法得到的地震热点内。门源M_S6.4地震前存在明显的地震热点，且地震热点相对集中，地震热点图像颜色由深变浅，震后消失；在祁连山北缘断裂附近，地震热点持续存在，分布较集中，地震热点图像颜色也呈由深变浅的过程。未来需关注祁连山北缘断裂，附近有发生强震的可能。     

Generalized Jacobi spectral Galerkin method for fractional pantograph differential equation

This work is concerned with the extension of the Jacobi spectral Galerkin method to a class of nonlinear fractional pantograph differential equations. First, the fractional differential equation is converted to a nonlinear Volterra integral equation with weakly singular kernel. Second, we analyze the existence and uniqueness of solutions for the obtained integral equation. Then, the Galerkin method is used for solving the equivalent integral equation. The error estimates for the proposed method are also investigated. Finally, illustrative examples are presented to confirm our theoretical analysis.     

铟镓共掺杂对n-ZnO纳米棒/p-GaN异质结生长行为和光电性能的影响

利用低温水热法在p-GaN薄膜上生长了铟(In)和镓(Ga)共掺杂的ZnO纳米棒。X射线衍射(XRD)、X射线光电子能谱(XPS)和X射线能量色谱仪(EDS)结果表明,In和Ga已固溶到ZnO晶格中。扫描电子显微镜(SEM)结果表明, ZnO纳米棒具有良好的c轴取向性,随着In和Ga共掺杂浓度的增加,纳米棒的直径减小,密度增加。XRD结果表明,In和Ga共掺杂引起ZnO晶格常数增大,导致(002)衍射峰向低角度方向偏移。同时,ZnO的光学性质受到In和Ga共掺杂的影响。与纯ZnO相比, 共掺杂ZnO纳米棒的紫外发射峰都出现轻微红移,这是表面共振和带隙重整效应综合作用的结果。I-V特性曲线表明,随着In和Ga共掺杂浓度的增加,n-ZnO纳米棒/p-GaN异质结具有更好的导电性。     

Chiral one‐dimensional hydrogen‐bonded architectures constructed from single‐enantiomer phosphoric triamides

The two single‐enantiomer phosphoric triamides N‐(2,6‐difluorobenzoyl)‐N′,N′′‐bis[(S)‐(−)‐α‐methylbenzyl]phosphoric triamide, [2,6‐F₂‐C₆H₃C(O)NH][(S)‐(−)‐(C₆H₅)CH(CH₃)NH]₂P(O), denoted L‐1 , and N‐(2,6‐difluorobenzoyl)‐N′,N′′‐bis[(R)‐(+)‐α‐methylbenzyl]phosphoric triamide, [2,6‐F₂‐C₆H₃C(O)NH][(R)‐(+)‐(C₆H₅)CH(CH₃)NH]₂P(O), denoted D‐1 , both C₂₃H₂₄F₂N₃O₂P, have been investigated. In their structures, chiral one‐dimensional hydrogen‐bonded architectures are formed along [100], mediated by relatively strong N—H…O(P) and N—H…O(C) hydrogen bonds. Both assemblies include the noncentrosymmetric graph‐set motifs R₂²(10), R₂¹(6) and C₂²(8), and the compounds crystallize in the chiral space group P1. Due to the data collection of L‐1 at 120 K and of D‐1 at 95 K, the unit‐cell dimensions and volume show a slight difference; the contraction in the volume of D‐1 with respect to that in L‐1 is about 0.3%. The asymmetric units of both structures consist of two independent phosphoric triamide molecules, with the main difference being seen in one of the torsion angles in the OPNHCH(CH₃)(C₆H₅) part. The Hirshfeld surface maps of these levo and dextro isomers are very similar; however, they are near mirror images of each other. For both structures, the full fingerprint plot of each symmetry‐independent molecule shows an almost asymmetric shape as a result of its different environment in the crystal packing. It is notable that NMR spectroscopy could distinguish between compounds L‐1 and D‐1 that have different relative stereocentres; however, the differences in chemical shifts between them were found to be about 0.02 to 0.001 ppm under calibrated temperature conditions. In each molecule, the two chiral parts are also different in NMR media, in which chemical shifts and P–H and P–C couplings have been studied.